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LONG-TERM  GOALS 

The  long-term  goals  of  this  research  have  been  to  quantify  the  incidence  of  surface  wave  breaking  and 
its  role  in  dissipating  wave  energy,  which  is  then  available  to  generate  currents  and  turbulence  in  the 
upper  ocean. 

OBJECTIVES 

The  objectives  of  this  project  were  to  develop  and  use  remote  and  in  situ  techniques  to  measure  the 
kinematics  and  dynamics  of  breaking  and  relate  them  to  the  dynamics  of  the  wave  field  and  the  surface 
mixed  layer.  Of  particular  interest  has  been  the  dissipation  of  surface  wave  energy  by  breaking. 

APPROACH 

The  approach  revolved  around  the  use  of  aerial  imaging  from  light  aircraft  in  the  coastal  zone  and 
quantitative  image  sequence  analysis  using  techniques  of  Particle  Imaging  Velocimetry  (PIV)  to 
measure  the  kinematics  of  whitecaps.  Kinematic  measurements  were  used  to  investigate  aspects  of 
surface  mixing  including  estimates  of  the  rate  of  mixing  of  the  surface  layer.  In  conjunction  with 
simple  scaling  arguments  to  relate  the  kinematics  of  breaking  to  the  dynamics,  and  laboratory 
experiments  on  breaking,  the  remote  data  were  used  to  investigate  wave  dissipation  and  the  momentum 
flux  from  waves  to  currents.  In  situ  measurements  of  waves  and  currents  to  compare  with  the  remote 
measurements  remains  a  topic  of  continuing  research. 

WORK  COMPLETED 

This  project  has  focused  on  the  Shoaling  Waves  Experiment  (SHOWEX)  conducted  off  Duck,  N. 
Carolina  in  the  fall  of  1999  with  other  supporting  measurements  off  Oahu,  Hawaii,  and  California. 
During  SHOWEX  our  imaging  system,  the  Modular  Aerial  Imaging  System  (MASS)  was  flown  on  the 
Long-EZ  by  Tim  Crawford  (deceased)  of  NOAA.  Using  the  data  acquired  from  that  experiment  we 
developed  algorithms  based  on  both  PIV  and  simpler  techniques  to  analyze  the  image  sequences  to 
measure  A(c)  dc,  the  mean  length  of  breaking  fronts  per  unit  area  of  ocean  traveling  with  speeds  in  the 
range  (c,  c+dc)  (Phillips,  1985).  This  statistic  forms  the  basis  of  a  hierarchy  of  moments  that  can  be 


used  to  investigate  the  kinematics  and  dynamics  of  breaking.  Figure  1  shows  examples  of  the  image 
data,  the  traditional  whitecap  coverage  and  an  example  of  the  PIV  processing  on  a  single  whitecap. 

Two  moorings  were  deployed  to  measure  near-surface  turbulence  (ADV  and  Dopbeam  sensors)  and 
temperature  fields  and  bubbles  entrained  by  breaking  waves.  Two  bottom  ADCP/pressure  packages 
were  deployed  to  measure  waves  and  currents  on  the  shelf.  One  mooring  suffered  a  failure  during  a 
hurricane  and  was  redeployed.  One  of  the  moorings  was  inexplicably  destroyed  by  the  Canadian  Coast 
Guard  vessel  Creed,  a  participant  in  the  experiment,  during  good  weather  with  unlimited  visibility. 

RESULTS 

MASS  data  from  SHOWEX  has  been  analyzed,  written  up  and  published  by  Nature  (Melville  & 
Matusov,  cover  article,  May  2, 2002).  We  have  shown  that  during  SHOWEX  A(c)dc  increases  like  the 
wind  speed  cubed,  Uio3,  and  decays  exponentially  with  the  breaker  speed,  c.  Since  the  breaking 
statistics  depend  on  the  directional  wave  spectrum,  and  since  there  is  no  universal  wave  spectrum  over 
the  full  range  of  breaking  waves,  we  do  not  expect  any  universal  shape  for  A(c)dc.  Simple  physical 
and  dimensional  arguments  provide  interpretations  of  the  various  moments  of  A(c).  For  example,  the 
first  moment  is  proportional  to  the  rate  of  mixing  of  the  surface  area  by  breaking,  the  third  is 
proportional  to  the  volume  mixed,  the  fourth  is  proportional  to  the  momentum  flux  from  waves  to 
currents  and  the  fifth  is  proportional  to  the  wave  dissipation.  These  moments  are  shown  in  Figure  2. 
Note  that  as  you  move  to  higher  moments  the  mode  of  the  distribution  moves  to  larger  values  of  c. 
Since  c  scales  approximately  with  the  phase  speed  of  the  breaking  waves,  this  implies  that  while  small- 
scale  waves  contribute  most  to  the  length  of  the  breaking  fronts,  longer  waves  are  responsible  for 
momentum  transfer  and  dissipation.  Since  the  vector  c  and  not  just  the  speed  are  measured,  we  can 
examine  the  directional  distribution  of  A(c)  and  its  moments  can  be  measured.  Examples  in  Figure  3 
show  that  for  the  SHOWEX  data,  A(c)  and  its  moments  are  approximately  symmetric  about  the  wind 
direction. 

The  in  situ  data  from  the  surface  moorings  shows  that  the  dissipation  of  turbulent  kinetic  energy  (TKE) 
measured  by  the  Dopbeam  correlates  with  the  modulation  of  the  surface  waves  or  wave  groups.  An 
example  of  this  data  is  shown  in  Figure  4.  Analysis  of  the  Dopbeam,  ADV  and  temperature  data 
continues. 

IMPACT/APPLICATIONS 

During  this  project  we  have  shown  that  methods  of  PIV  developed  in  the  laboratory  can  be  used  to 
analyze  airborne  imagery  to  measure  the  kinematics  and,  with  supporting  assumptions,  the  dynamics  of 
wave  breaking.  These  results  are  novel  and  are  the  consequence  of  the  convergence  of  technical 
developments  in  the  areas  of  imaging,  position/motion  detection  (DGPS  and  inertial  systems)  and 
imaging  software.  While  we  have  concentrated  on  visible  imagery,  the  same  approach  could  be  taken 
with  IR  and  hyperspectral  imagery  to  study  the  kinematics  of  different  fields  including  temperature. 

The  in  situ  data  shows  great  promise  for  making  acoustical  measurements  of  ocean  microstructure 
directly  in  the  wavenumber  domain  (Veron  &  Melville,  1999).  This  avoids  the  use  of  Taylor’s 
hypothesis  to  move  from  the  frequency  domain  to  wavenumber  domain;  an  approach  that  has  been 
fraught  with  problems  in  the  surface  wave  zone  which  is  dominated  by  the  orbital  motion  of  the  waves. 


TRANSITIONS 


The  imaging  system  has  been  transitioned  to  the  measurement  of  breaking  and  spray  generation  in 
hurricanes  with  support  from  ONR’s  CBLAST  program. 

RELATED  PROJECTS 

Related  projects  are  described  on  our  web  page  at  http://airsea.sio.ucsd.edu. 
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Distance  off  shore,  m 


F  IG.Ia-c:  Images  showing  approximately  160  m  x  160m  of  the  sea  surface  for  mean  wind 
speeds,  U-io,  of  7.2,  9.8  and  13.6  m  s'1,  respectively,  in  SHOWEX.  Note  the  increasing 
density  of  whitecaps  with  Uio  and  their  random  shapes,  d  Whitecap  coverage  as  the  aircraft 
flew  offshore.  Note  the  large  fluctuations  at  spatial  scales  of  O  (1- 10)  km  that  may  be  due  to 
both  local  wind  and  wave  modulations,  e:  The  result  of  processing  a  single  whitecap  using 
the  PIV  technique  in  which  the  normal  velocity  of  the  boundary  of  the  whitecap  is  resolved. 
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FIG.  2a,  Binned  measurements  of  A(c)  weighted  by  Ui0'3  showing  that  the  average  length  of 
breaking  fronts  in  (c,  c+dc)  per  unit  area  of  sea  surface  increases  like  U-io3  and  decays 
exponentially  with  c:  A(c)  =  S.SxIO^.e'0  640. 

b.  The  weighted  first  moment  of  A(c)  which  corresponds  to  the  fractional  area  swept  out 
(turned  over)  by  breakers  in  the  speed  range  (c,c+dc)  per  unit  time. 

c.  The  weighted  fourth  moment  of  A(c),  c4A(c)w,  which  corresponds  to  the  momentum  flux 
from  waves  to  currents  due  to  whitecaps,  with  the  subscript  “w”  referring  to  the  component  in 
the  wind  direction. 

d.  The  weighted  fifth  moment  of  A(c),  which  corresponds  to  the  energy  lost  from  the  wave 
field  due  to  breaking:  “wave  dissipation". 
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Figure  3a.  Distributions  of  A(c),  with  respect  to  the  northerly  and  easterly  components  of  c 
(cn.Ce),  along  with  the  average  wind  direction  (arrow).  Note  that  the  dominant  orientation  is 
close  to  the  wind  direction. 

b.  The  corresponding  momentum  flux  from  waves  to  currents  due  to  breaking  for  the  three 
wind-speeds.  Note  that  as  the  wind  increases  the  region  of  significant  momentum  flux 
becomes  approximately  symmetrical  about  the  downwind  direction  (c.f.  Fig.  2b). 

c.  The  corresponding  distribution  of  the  dissipation  with  (cN,cE). 


FIG.4.  Top:  Inertial  estimates  of  dissipation  from  the  Dopbeam  acoustic  Doppler  near  the 
surface  in  SHOWEX. 

Bottom:  Simultaneous  measurements  of  the  vertical  component  of  orbital  velocity  measured 
with  an  ADV  on  the  same  mooring.  Note  the  corelation  between  the  modulation  of  the  wave 
field  and  the  dissipation,  especially  around  45-50  s. 


